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The ttN s-wave scattering lengths have been inferred from a joint analysis of the pionic hydrogen 
and the pionic deuterium x-ray data using a non-relativistic approach in which the ttN interaction 
is simulated by a short-ranged potential. This potential is assumed to be isospin invariant and its 
range, the same for isospin 1=3/2 and 1=1/2, is regarded as a free parameter. The proposed model 
admits an exact solution of the pionic hydrogen bound state problem, i.e. the ttN scattering lengths 
can be expressed analytically in terms of the range parameter and the shift (e) and width (F) of 
the Is level of the pionic hydrogen. We demonstrate that for small shifts and short ranges from the 
exact expression one retrieves the standard range independent Deser-Trueman formula. The 7rd 
scattering length has been calculated exactly by solving the Faddeev equations and also by using a 
static approximation. It has been shown that the same very accurate static formula for vrd scattering 
length can be derived (i) from a set of boundary conditions; (ii) by a reduction of Faddeev equations; 
and (iii) through a summation of Feynman diagrams. By imposing the requirement that the 7rd 
scattering length, resulting from Faddeev-type calculation, be in agreement with pionic deuterium 
data, we obtain bounds on the ttN scattering lengths. The dominant source of uncertainty on the 
deduced values of the ttN scattering lengths are the experimental errors in the pionic hydrogen 
data. 
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I. INTRODUCTION 



The determination of low-energy pion-nucleon (ttN) parameters has been the focus of much theoretical and ex- 
perimental efforts. The s-wave ttN scattering lengths are of particular importance serving as testing ground for 
various theoretical considerations. In addition to that, their isovector combination provides input in the Goldberger- 
Miyazawa-Oehme Q sum rule to be used to extract the ttNN coupling constant. In recent years major advances have 
been made in the experimental and theoretical investigation of the ttN system. With the advent of meson factories 
(LAMPF, PSI and TRIUMF) and the corresponding influx of the new high accuracy ttN scattering data considerable 
progress has been achieved in the ttN phase shift analyses ||-|| providing means to examine even such subtleties 
as isospin symmetry breaking effects [^,^,^ Recently, the ttN scattering experiments have been complemented by 
high quality pionic x-ray measurements performed, both on pionic hydrogen and on pionic deuterium The 
measurements of the shifts and widths in the Is levels in these atomic systems, resulting from strong ttN interac- 
tion, allows to extract directly the corresponding scattering lengths, i.e. a^^p and a^^d, respectively. Therefore, the 
new x-ray data constitute an independent source of information on the low-energy ttN scattering parameters. On 
the theoretical side, the physical quantities bearing on the low-energy ttN interaction have now become accessible 
to calculations conducted within quantum chromodynamics (QCD). Since QCD is known to be highly non- 
perturbativc at low energies, its low-energy implementation has been based instead on a chiral perturbation theory 
in which the effective Lagrangian is expanded in increasing powers of derivatives in meson fields and quark masses. 
This approach in practice involves a Taylor expansion in the meson four-momenta and therefore it may be expected 
that the lower the energy is, the more accurate are the predictions. In this context, the precise knowledge of the 
experimental values of the low-energy ttN scattering parameters is essential for further development of the theory. 

The purpose of this work is to extract the s-wave ttN scattering lengths using exclusively the pionic hydrogen 
and pionic deuterium x-ray data. The key reason for proceeding along this route is that the low-energy regime can 
be thereby investigated without recourse to scattering data and there is no danger that the low-energy parameters 
have been largely determined by the data at high energies. Our treatment is purely phenomenological based on an 
isospin invariant potential model and we wish to clarify at the onset that this approach relinquishes any pretense of 
being a theory in favour of practicable calculational scheme. The investigation has two parts. In part one we take 
as our input the values of the ttN scattering lengths determined previously from pionic hydrogen data and use them 
in a microscopic calculation of the 7rd scattering length. The latter has not been measured directly in a scattering 
experiment but may be extracted from the pionic deuterium x-ray data by applying the Deser-Trueman formula 
|ll| . It is an empirical fact that the ttN scattering lengths are small as compared with the deuteron size and it 
has been a common practice JT2[ to use the multiple scattering expansion for calculating the 7rd scattering length. 
Since this series rapidly converges, what has been confirmed by early Faddeev calculations |13-15|, in the past with 
the poorly known 7rd scattering length there was little incentive to go beyond the second order (for a review, cf. 
[fi6|-|l8[). At present, the experimental error on 7rd scattering length is at the level of 2% and the adequacy of the 



second order formula might be questionated. Strictly speaking, a truncation of the multiple scattering series can 
really only receive its justification when we actually quantify the magnitude of the higher-order terms to establish 
whether they are truly negligible. This question is examined in detail in this paper and the Trd zero-energy scattering 
problem is solved exactly within a three-body formalism by introducing a zero-range model to simulate the ttN s-wave 
interaction. One advantageous feature of this model is that it allows to obtain an analytic solution of the three-body 
problem in the static approximation. We demonstrate that the static solution can be obtained either by reduction of 
the Faddeev equations, or by imposing a suitable set of boundary conditions, or finally by performing a summation 
of Feynman diagrams. All three methods converge to the same analytic formula expressing the 7rd scattering length 
in terms of the ttN scattering lengths. Static solution in coordinate space is very appealing and helps to develop an 
intuitive picture of how the individual ttN amplitudes contribute to build up the ird scattering length. By solving 
numerically the Faddeev equations we show that the accuracy of the static approximation is comparable with the 
present experimental uncertainty on a^d- In order to find out what the pionic deuterium data can teach us about 
the ttN scattering lengths, the 7rd scattering lengths obtained as a solution of the Faddeev equations is compared 
with experiment. It turns out that the three-body calculation is in agreement with experiment only when the input 
ttN scattering lengths belong to a relatively small subset of values that are consistent with pionic hydrogen data. 
The ttN scattering lengths that belong to this subset simultaneously satisfy the constraints imposed by the pionic 
hydrogen and pionic deuterium data. 

In part two of the present work we introduce explicitly a range parameter in order to examine the validity of 
the zero-range model. To achieve this goal it is essential to devise a simple and transparent representation of the 
ttN interaction in which the two-body scattering problem with and without Coulomb interaction admits an analytic 
solution and we show that a two-channel isospin invariant separable potential lends itself to that end. Moreover, 
within this representation the exact bound state condition appropriate for the pionic hydrogen problem takes also an 
analytic form. The latter being a single complex constraint, is equivalent to two real equations that can be explicitly 
solved and as a result the ttN scattering lengths are obtained as functions of the range parameter together with 
the Is level shift and width in the pionic hydrogen. In particular, when the level shift is small as compared with 
the Coulomb energy and the range of the interaction is small in comparison with the Bohr radius, from the exact 
bound state condition we retrieve the Deser-Trueman formula (independent of the range parameter). Regarding the 
range as a free parameter we are able to extend the zero-range model and by varying this parameter in physically 
reasonable limits we find the results to be insensitive to the value of the range. The uncertainty on the ttN scattering 
length caused by the lack of knowledge of the range is much smaller than that resulting from the experimental errors 
on the pionic hydrogen level shift and width. 

The organization of this paper is as follows. In Sec. || we develop a zero-range model and review various derivations 
leading to the static solution of the 7rd scattering problem. The accuracy of the static solution is examined by 
comparing it with the solution of the Faddeev equations. We infer isoscalar and isovector ttN scattering lengths that 



are consistent with both, pionic hydrogen and pionic deuterium data. In Sec. [II we lift the zero-range limitation by 



introducing a finite range into our formalism. We present an exact treatment of the pionic hydrogen and we derive 
Deser-Trueman formula for that particular case. The 7rd scattering length obtained from the solution of the Faddeev 
equation is compared with experiment. Finally, the results are summarized in Sec. 



II. ZERO-RANGE MODEL 



The central issue we wish to address in this section is how to construct a theoretical framework in which we can 
use the pionic deuterium data to gain information on the ttN scattering lengths. The measurement of the shift and 
the width of the Is level in pionic deuterium presents us with the value of Trd scattering length ajrd. The latter 
quantity is defined as the elastic Trd scattering amplitude evaluated at zero kinetic energy of the incident pion. This 
amplitude is necessarily complex because absorption reaction channels are open even at the very threshold. The most 
important of them is the Tr~c? —f nn reaction, and to a lesser extend the radiative absorption t:~ d — > ^nn channel. 
In principle, there would be also the charge- exchange break-up channel Tr~d — > n'^nn that is open at threshold but 
this process is strongly suppressed by the centrifugal barrier. Indeed, with s-wave ttN interaction there is no spin-flip 
possible so that for the two neutrons the ^Sq state is not available, whereas the ■^S'l state is forbidden and they have 
to be produced in higher partial waves. On the whole, however, the absorptive effects are not large at threshold, 
judging from the magnitude of the imaginary part of the Trd scattering length which empirically constitutes only 
about a quarter of the real part of a^^d- Strictly speaking, the absorptive processes contribute to both, the real and 
the imaginary part of ajrci but in the following we are going to ignore the absorptive corrections to the real part of 
a^d- Disregarding the absorptive processes, we shall concentrate our attention on a microscopic calculation of a,rd 
and in order to be able to solve the ensuing three-body problem we introduce a potential description of the TrN 
interaction to be used in the appropriate Faddeev equations. 

In order to facilitate the discussion of the Faddeev approach, it is instructive to take the static model as our point 
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of departure. The attractive feature of the static model is that it is much easier to develop and to compute since 
the final result for pion-deuteron scattering length takes the form of a single analytic formula that does not require 
off-shell information. Moreover, in our case the latter model also happens to be extremely good approximation to the 
full solution of the three-body problem. The earliest version of a static model, due to Brueckner was based on 
the fixed scatterer concept and ignored all isospin complications. Here, we wish to make it somewhat more realistic 
introducing as our dynamical framework a set of appropriate boundary conditions, but on the other hand, we are 
prepared to content ourselves with a theory that has isospin invariant point like interactions. Labeling the pion as 
1 and the nucleons as 2 and 3, the boundary conditions representing the zero-range tt-N interaction taking place on 
nucleon i where i = 2,3, may be written as 



lim \xi~Xi\'^{xi,X2,X3) ^ {^j,/77i){bo + biI -T^) lim -— - a;i|^'(a;i, a;2, 0:3) (1) 

where the bar denotes an average over directions xi — Xi what is equivalent to projecting out the s-wave component 
of the wave function 5*, and the boundary condition is to be imposed for each of the two nucleons. The vectors 
I and T are, respectively, the pion and the nucleon isospin operators, whereas &o and hi denote the isoscalar and 
isovector tt-N scattering lengths, is the tt-N reduced mass and m is the pion mass. In the following we choose the 
cm. of the two nucleons as the origin of the coordinate system, i.e. we set X2 — ■^r and x^ = "i'" with r being the 
nucleon-nucleon separation vector. The pion vector in this Jacobi coordinate system will be denoted as p. When 
the wave function \E'(r,p) describing the ttNN system for the case of 7t~ scattered off the deuteron is known, the 
amplitude leading to the final state with asymptotic wave function is — ($/|T^|\E') where V denotes the potentials 
that have been taken out in the derivation of <&/. For elastic scattering "I>/(p, r) = exp (zp' • p)ipd{r) where p' is 
the momentum of the outgoing pion, V-'d is the deuteron wave function and V is the sum of the two ttN potentials 
as asymptotically there is no 7r-deuteron interaction. Although in our formalism we never needed ttN potentials and 
the ttN interaction is represented by the boundary condition (|^), it is in fact possible to give a formal expression for 
such potential (cf. [Q) and for the operator V we take 

V^ip, r) = |(6o + &i J • T2) d{p - ir) A|p _ i^-] + (bo + hI- T3) S{p + ir) A|p + r). (2) 

Denoting the incident pion momentum as p and making use of the boundary conditions (0) in (|^), the 7r-d elastic 
scattering amplitude f{p\p) takes the form 

fip',P)^^ J e-'P'P^^M {6ip ~ ^r)\p - ^r\ + 6{p + ^r)\p + ^r\} ^ip,r)d^p d\ (3) 

where v is n-d reduced mass. Given the elastic n-d scattering amplitude (^), the n-d scattering length follows 
immediately from 

a.d-/(0,0). (4) 

With the TT-N interaction assumed to be isospin invariant, it will be convenient for us to adopt an isospin notation. 
For the initial 7r~-d system, the isotopic spin wave function has the form 

Xa = 7r"^(p2"-3 - "2^3), (5) 

where the symbols p, n, it~ in (|^) stand for the isospin wave functions of the corresponding particles. The wave 
function is antisymmetric in the nucleon labels, as appropriate for the state where the isospin of the two-nucleon 
subsystem equals zero. As a result of the interaction, the two nucleons can undergo a transition to a symmetric 
configuration corresponding to /23 — 1 and we shall need also a function that is symmetric under two-nucleon 
permutation 

Xs = ^TT" (p2n-3 + n2P3) - ^7r° 722713. (6) 

Since our interest here is confined to s-wave interactions, no spin flip is possible and therefore the spin part of the 
wave function does not change. Regarding the nucleons as fixed scattering centers, we may anticipate that the wave 
function '9(p,r) for the full system of the target nucleons and the meson will take the approximate form 



^(p, r) = e'P P Ud{r) Xa + A(r) 



exp {ip\p- \r\) ^ cxp(zp|p+ ^r|) 
\p-\r\ \p+\r\ 



+ X{r) 



exp {ip\p — \r\) exp {ip\p 



2' \i _ < 2' 

ir| |p+ir| 



Xa + 

(7) 
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where Ud is the spatial part of deuteron wave function that includes also the deuteron spin and in particular may 
contain also the D-component. The projectile enters with momentum p and in the initial asymptotic region the pion 
and the target have separate wave functions (a plane wave and Wd(r), respectively) and the propagation from one 
scattering center to another is described by a superposition of spherical waves. The hitherto unknown amplitudes 
denoted in respectively, as A{r) and X(r) multiplying these outgoing waves emitted by the two centers account 
for the multiple scattering phenomena. They will be determined from the boundary conditions (|l|). To satisfy Pauli 
principle the wave function (j^) must be antisymmetric in the two nucleon variables. This implies that we have to 
stipulate that the coefhcients A{r) and X{r) are even under permutation of the nucleons, i.e. they must be invariant 
under the reflections r — > — r. For zero-energy scattering considered in this work, however, this is always the case 
because A{r) and X{r) depend then only upon the magnitude of r. It is worth noting that the wave function 
includes explicitly virtual charge exchange amplitude X{r). Since our interest here is confined to zero-energy 
scattering, in the following we take p = in (|^). Equations for the functions A{r) and X{r) may be obtained by 
substituting (|^) in (|l|) for i — 2 and equating the coefficients multiplying the same isospin functions. With two 
different isospin functions we obtain two equations and this procedure determines uniquely A(r) and X(r). Owing 
to the proper antisymmetrization of our wave function the boundary condition for i = 3 will be then automatically 
satisfied. The equations obtained from are 

A{r) = boudir) + (bo/r) A{r) + V2{h/r) X{r), (8a) 
-X{r) = V2biUdir) + V2{bi/r) A(r) + (bo + h)/rX{r). (8b) 

In (^) we introduced the abbreviation bj = (1 + m/M)bj where M is the nucleon mass. The n-d scattering length is 
given by the overlap integral 

a^d - (2i^/m) J Mrf(r)U(r)d3r (9) 

where A{r) is the solution of (^) 

Air) = bo + ibo + b,)ibo-2b,)/r 

l-foiA--(feo + &i)(&o-26i)/r2 

Using (|o|) in and expanding A{r) in powers of the ttN scattering lengths, we retrieve the well known second 
order formula for the 7r-d scattering length (cf. p8[) 



(2) _ 2zy 



(6g-2fo?)(i 



(11) 

where the expectation value is taken with respect to the deuteron wave function. As advertised at the beginning of 



this section, formula (10) provides a complete solution of the problem. To examine the accuracy of the static formula 
we have to compare it with the exact solution of the three-body problem. The latter will be obtained by solving the 
Faddeev equations on which we now embark. 

To solve the Faddeev equations it will be convenient for us to work in momentum space. Introducing the Faddeev 
partitions, we write the three-body wave function as 

* = ^(1) {qi , fci) + (g^^ ^ ^(3) (g^^ fc^)^ (12) 

where qi denotes the relative momentum of the (23) pair whereas fci is the cm. momentum of particle 1 and cyclic 
permutations are implied. To obtain Faddeev equations for the amplitudes, the different partitions are written as 
(cf. lH) 

^^'\q, k) - (27r)30(q)<5(fc - p)xa + [F{q, k)xa + G{q, k)xs]/{E - q'/M - ^ /2v)- (13a) 
xl;^^\q, k) = [A{-q, k)xa - X{-q, k)xs]/{E - /2n - ^12^^); (13b) 
V'(3)(q, k) = [A(q, k)xa + Xiq, k)xs]/{E - 12^ - l2vt,Y (13c) 

where i/jv is the reduced mass of the nucleon and that of the ttN pair, E is the cm. three particle kinetic energy and 
(/>(q') is the deuteron wave function in the momentum space. In (|l^) we have introduced four scattering amplitudes 
F{q,k),G{q,k), A{q,k) and X{q,k). However, the amplitude G{q,k) to be non-zero requires at least p-wave NN 
interaction and therefore will be excluded from our considerations, while the three remaining amplitudes will be 
determined from the Faddeev equations. It is evident from ( p^ ) that under the P23 permutation ■f/;^^-' — ■(/j^^) and 
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^(2) ^ -tp'^^\ SO that the total wave function is, as required, antisymmetric in the nucleon labels. Assuming exact 
isospin conservation, we can write the Faddeev equations 



F{q,k) 



d3 fc' {q\tiE - ei2v)\\k + k') + {q\t(E - e/2v)\ 



(27r)3 



E~{k + ^ik'/MYl2^i - fc'2/2iyjv 



(14a) 



A{q,k) 



to{E- 



2vN 



M 



k+p) 0(fc+ ip) + 



d^k' {q\to{E - ky2iyN)\f^k/M + k') 



(27r)3 E-{k + \k'f/M - k''^/2v 



^k')^/M-k'^ 
d^k' {q\ta{E-k'^/2vN)\~ iik/m-k') 



F{-k - ife',fc') 



(27r)3 E- 



{k + fik'/my/2fi- k'^/2i^N ^ 



/- f d3fc^ {q\h{E-e/2yN)\~^^k/m~k') _ a, 

J (2^)3 E-{k + iJLk'/mY/2^i-k'^/2vN ^ m ' ' 



(14b) 



-X{k,k) 



V2(k 



■^k 



MH- {kME-ky2^^)\f^k/M + k') , 
+'' 7(2^)3 - (fc + ifc') VAf - fc'V2^. ^ 



) + 



d3 A:' (fc|[to(£^ - k^/2i^N) ~h{E- ^ I2vn) 



(2^)3 

, /2 f d^ (k\ti{E_Pl 
J (27r)3 E-{k + fik' 



(fc + nk'/m)^ /2fi - fc'2 /2!yAr 
/2i^jv)| — fik/m — fc') 



^^/"l- ^x(-fc-iifc',fc') 

m 



/mY/2^i-k'^/2vN ^ m^''^'^' 



(14c) 



where in (|T^) (q'|<(_E)|<7) is the NN scattering t-matrix for zero isospin and {q'\tj{E)\q) are, respectively, the isoscalar 
{i — 0) and isovector (j = 1) ttN scattering t-matrices. The elastic scattering amplitude is given by the expression 



/(P',P) = lim 



47r 



F{q.p') 



d\ 

E-q^lM-p'^l2v (27r)3 



(15) 



and the scattering length is obtained from (Q). We can use ( |14a| ) to eliminate F{q,k) in ( p^ in favour of the 
amplitude A{q,k). In the NN scattering matrices occurring in (p_4D, as a result of the limiting procedure, only the 
deuteron pole contributes and scattering length is given as an overlap integral 



(16) 



The above formula is analogous to (||), and, in fact, the static approximation results (||)-(|l^) could have been derived 
from the Faddeev formalism. In order to demonstrate that (^)-(|l^) follow from ( p^ ) we note that when the nucleons 
are static they are not supposed to scatter (t 0) and the amplitude F{q, fc) drops out in (14b) and (14c) so that 
we are left with only two coupled integral equations. When the underlying forces are of zero range, the off-shell ttN 
scattering amplitudes can be simplified, and in that case 



-(Mm) bj/il + ^b, 



J =0,1; 



where = 2^B and B is the binding energy of the deuteron. The important consequence of the zero-range 
assumption, apparent from the above formula, is that the t-matrices become independent upon the off-shell momenta. 
Therefore, the amplitudes A{q, fc) and X{q, fc) will be functions of one variable only and it will be convenient for us 
to introduce a notation that emphasizes that fact, setting A{q,k) — --{m/2Tr)A{k) and X{q,k) — — (m/27r) A'(fc) , 
wher e A( k) and X{k) are two, hitherto unknown amplitudes. With static nucleons, the energy denominators in ( |l4b| ) 
and (14c) become all equal to ~B — (fc' + k)^ /2m and we end up with the following set of integral equations for the 
amphtudes A{k) and X{k): 



A{k) = bQ(j){k) + Att bo 



d^k' 



Aik') 



(27r)3 k2 + (fc' + fc) 



X{k') 



(fc' + fe)2 



(17a) 
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where 

^bj{l + m/M)/{l + Kbj); j = 0,l. (18) 
The above set of integral equations can be immediately solved by introducing the Fourier transform 

A{r) = / e''^M(/c)d3A: (19) 



together with a similar relationship for X(k) and (f>{k) and using the well known formula 



K2 + (fc + fc')2 

In order to solve ( p7| ) we multiply the latter equations by e**"" and subsequently integrate them over fc. As a result, 
we obtain a set of two algebraic equations for A(r) and X[t) that differ from (@) only by exp {—Kr)/r replacing l/r 
and bj replacing bj. Since (|l^ ) goes over into we are led to the extension of the static formula 



. bo + ibo + h)ibo-2h)e-^yr 

A(r) = Ud[r). (20) 

1 - 5ie-«7r - (6o + 6i)(5o - 26i)e-2«7r2 

This formula is to be used in (|^) but now accounts for the binding energy correction. 

Concluding our discussion of the static model we wish to recall that a closed form expression for vrd scattering 
length has been also obtained by effecting an explicit summation of Feynman diagrams and the most complete 
treatment can be found in Ref. |2l| . The ultimate static formula for a^d, that takes into account isospin degree of 
freedom, given in |2l] is rather complicated and at first sight appears to be different from (pO|). However, a closer 
inspection reveals that the authors of Ref. apparently did not realize that their fractional formula for a-^d could 
have been significantly simplified because a common factor equal 

1 + bie-'^Vr - (feo + &i)(6o - 2bi)e-^^'' /r'' 

may be pulled out both, from the numerator, as well as from the denominator and eventually drops out. Indeed, when 
the redundant factor has been cancelled, the resulting expression is identical with (|2^). Therefore, when binding 
corrections are disregarded, this approach reproduces the static model result ( p^ and it is reassuring that in this 
case all three methods give the same answer. 

To improve upon the static model one needs a numerical solution of the Faddeev equations and in the following, 
similarly as in the previous calculations |ll3|-|l5|, in order to reduce the computational effort, all the pairwisc interac- 
tions invoked will be represented by rank-one separable potentials. The ttN s-wave interaction is taken in the form of 
a standard Yamaguchi potential with the same form factor in both isospin states. Since the inverse range parameter 
/3 that enters that form factor is not known, similarly as before, we consider the zero-range limit, i.e. (3 oo. For 
an assigned value of /3, the strength parameter of the potential may be eliminated in favour of the scattering length 
and the appropriate s-wave t- matrices, are 

{k\t,iE)\k') = - — —^^ l-,,pb,{l-2ip/f3){l-w/P)-^ 1 + ^^^^ 



where p — y/2iiE and j = 0, 1 and it is evident from (^) that the zero-range limit can be effected. When the 
nucleon motion is taken into account, the p-wave ttN interaction gives contribution to the 7rd scattering amplitude 
even at threshold. Therefore, in addition to s-wave, we are going to include also the p-wave interaction, limiting 
ourselves only to the P33 wave as in that case both the strength and the statistical weight are dominant rendering 
the remaining p-waves negligible. The corresponding p-wave form factor of the form 

9A{k)^k/{k' + pl) 

has been adopted from jlij with /3a — 5.33 fm^^ where the depth of the separable potential can be adjusted to the 
experimentally known value of the P33 scattering volume taken to be 0.64 fm^. It is well known that with the above 
form, the shape of the delta resonance cannot be well reproduced but this is less important here, the essential thing is 
to have the P33 amplitude at threshold correctly reproduced. Besides, the p-wave constitutes only a small correction 
and using a more complicated model does not seem to be currently justified. For the NN interaction we use two 
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separable models: a simple Hulthen-Yamaguchi potential with inverse range parameter equal Pn ~ 6.01162 \/MB 
whose strength is fixed by the deuteron binding energy, and the PEST potential constructed in Ref. |22| with a more 
sophisticated form factor of the form 



where the parameters Ci and (3i have been tabulated in Ref. |22|. This potential has been devised in such a way 
that the corresponding NN half-off-shell T-matrix has the same behaviour as that of the Paris potential |^ . This 
separable replica of the Paris potential takes into account the short range repulsion that is absent in the Yamaguchi 
potential yet retaining the simplicity of the latter. 

Using standard partial wave projections the Faddeev equations (|l4) can be reduced to a system of four coupled 
inhomogeneous integral equations in a single variable that are amenab e for numerical treatment. In the actual prac- 
tice, in order to cross-check our numerical procedures, we used two independent methods of solving these equations. 
The direct method introduces an integration mesh what allows us to replace integrals by sums so that the integral 
equations take the form of a system of linear algebraic equations easily solvable by standard methods. The second 
method solves the system of integral equations by successive iterations. The iterative procedure is equivalent to 
a power expansion in ttN scattering lengths what allows tracing down the contribution from the different orders. 
Since the scattering lengths are rather small, as compared with the deuteron size, the iterative sequence proves to 
be rapidly convergent. 

The experimental 7r-d scattering length has been extracted form the Is level shift in pionic deuterium by using the 
Deser-Trueman |ll[ formula. Therefore, the extracted quantity is in fact the Coulomb corrected scattering length, 
denoted hereafter as a^^, and before confronting the calculated pion-deuteron scattering length with experiment 
one needs the experimental value of a^^, i-S. the purely nuclear scattering length. Of course. Coulomb correction 
could be anticipated to be very small but since the experimental errors are also small, it is of interest to give some 
quantitative estimate of the Coulomb correction. In principle, for calculating the latter one needs to know the 
pion-deuteron nuclear potential responsible for the level shift. This potential is not known but with the zero-range 
potential simulating the ttN interaction in the first approximation it is reasonable to expect that the effective potential 
is proportional to the nuclear density, so that the shape of the nuclear potential is given by the square of the deuteron 
wave function U(i(2r)^. Still, the depth is not known, but on the nuclear scale this potential must be rather weak 
because the experimental value of a^^ is quite small. Therefore, to quantify the value of the ratio of a^Ja^^d it is 
sufficient to keep only the first order terms in the nuclear potential. Since in this case the potential depth drops out, 
we are led to the formula 

a%, ^ /-.^(2r)^0o(O,r)^dr ^ 
a^rd /g°°Md(2r)2r2dr 

where 0^ (fc, r) denotes the regular Coulomb wave function that for zero-momentum (k=0) and zero orbital momentum 
(i?=0), simplifies to the form 



00 (0, r) = r Ji (2\j2va.r) j^Phnxr (24) 



where a is the fine structure constant and J\(x) denotes the Bessel function. Expanding ( ^31) in powers of a, we 
obtain quite adequate first order formula a^^/oTrd = 1 ~ oiv (r) where the expectation value is with respect to the 
deuteron wave function. We have checked that for a variety of deuteron wave functions the calculated ratio (^3|) has 
ben very stable and its numerical value is 0.985. Using this number together with the experimental value of a^^ 

= [-(2.61 ± 0.05) -I- I (0.63 ± 0.07)] x 10~V™^ 

taken from Ref. where m^r is the mass of the charged pion, we deduce the value of the purely nuclear 7rd scattering 
length 

a^d = (-2.65 ± 0.05) x IQ-^/m^, (25) 

and hereafter the above number will be referred to as the " experimental" 7rd scattering length in which all absorptive 
effects have been neglected. 

Adopting the zero-range model of the ttN interaction, for calculating the Trd scattering length one needs as input 
just the isoscalar and the isovector ttN scattering scattering lengths. The values of bo and bi that have been extracted 
from the pionic hydrogen data in Ref. ^ , are 

&o = -(0.22 ±0.43) X IO^V^^tt; 6i = -(9.05 ± 0.42) x lO^V^^r, (26) 
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where the quoted uncertainty comprises the experimental errors together with the uncertainty introduced by applying 
a specific procedure that allows to deduce and bi from the measured x-ray spectra. The theoretical uncertainty 
is quoted to be about twice as large as the experimental error. Besides, the errors on bo and on bi are strongly 
correlated. 

Using (|2^) as our input, we have calculated the ird scattering length and the results are presented in Table 
H All entries are doubled because we employ two models of NN interaction: the numbers without brackets have 
been obtained using PEST wave function and, respectively, the bracketed quantities correspond to the Yamaguchi 
potential. For each set of input values of (60, &i) we computed a^^ using five different methods discussed before, 
beginning from the simplest second order formula ( pT| ) , through the static model (|l^) and (|2^) , up to the full Faddeev 
calculation without, and, with A, respectively. The results of the Faddeev calculation with s-wave interaction only 
(without A) constitute a benchmark for the various approximations. Contrary to what has been often claimed in 
the literature, the second order formula is insufficient as the error incurred is roughly four times bigger than the 
present experimental uncertainty on a^^. It is apparent from Table | that the closest to Faddeev result is in all cases 
the static model (|lO[). The accuracy of the latter is very good, the error being always below 2%. By contrast, the 
performance of the implementation (^) of static model is rather disappointing, especially that from formula ( po| ) 
containing binding energy correction, one might expect further improvement. Nevertheless, the numbers show just 
the opposite, that in fact the included corrections go in the wrong direction worsening the results so much that even 
the second order formula proves to be more accurate. Of course, it is not just the binding energy correction that is 
responsible for the difference between the static model and the Faddeev result, as only the latter properly accounts 
for the nucleon recoil. However, the lions share of the recoil correction seems to be cancelled with the binding energy 
correction and this cancellation explains the success of the static formula (|l^) containing neither of these corrections. 
An explicit demonstration that, at least to the second order, such mechanism is at work can be found in Ref. p^ . 

Since the static model (|o|) proves to be so accurate for Yamaguchi and PEST models of the NN interaction, we 
took advantage of this fact, using it to examine more realistic NN potentials containing also the D-wave part. The 
results of our computations are displayed in table || where we compare the two separable models (Hulthen- Yamaguchi 
and PEST), used in Faddeev calculations, with two popular local potentials (Paris and Bonn |^^). As expected, 
the PEST wave function results are indeed very close to those obtained with Paris wave function despite the lack of 
the D-component in the PEST wave function. Therefore, neglecting the D-wave in the Faddeev calculation does not 
appear to be a serious omission. It is also gratifying that PEST, Paris and Bonn models give very similar results. 

In table III we present the values of 7rd scattering length obtained in result of iterative solution of the Faddeev 
equations. Since for zero-range ttN interaction there is no additional suppression due to the ttN form factor, the 
rate of convergence is somewhat slower but the converged result is obtained in less than 10 iterations. We give 
a^d values calculated with and without p-wave ttN interaction what allows to evaluate the p-wave contribution in 
each order. For Yamaguchi NN interaction the p-wave correction in the first order is quite large and contributes 
0.47 X lO^^/m^. The p-wave contribution to the second order (called sp-term in Ref. [Q) has opposite sign and 
equals —0.35 x lO^^/m^. In general, the net effect of the p-wave interaction on the converged result is reduced 
owing to the destructive interference between repulsive s-waves and attractive p-waves, amounting in total only 
0.29 X 10~^/m7r. Similar features are observed for the PEST model but since the convergence rate is faster, the 
higher order corrections are suppressed and the interference effects seem to be smaller, i.e. the first order p-wave 
correction is 0.45 x 10~'^ /m-^ while the corresponding correction to the converged result is 0.39 x 10~^/m7r. 

It is apparent form table ^ that the calculated 7rd scattering length values are rather sensitive to the input values 
of (60,61) and therefore it is not so easy to see when the calculation agrees with experiment. To facilitate the 
comparison with experiment the values of a^^ resulting from Faddeev calculation (PEST with A) and displayed in 
table II have been represented analytically using bilinear interpolation on a grid in the (60,61) plane. Then, given 
the interpolating polynomial, we equated it to the experimental value of a.„d, adding or subtracting the experimental 
error. This procedure gave us two constraints of algebraic form in the (60, 61) variables, readily solvable with respect 
to one of these variables. The two functions obtained this way may be plotted in the (60, 61) plane where, as shown 
in Fig. they set the boundary of the tilted band representing the one standard deviation constraint imposed by 
the nd scattering length deduced form pionic deuterium data. The rectangle in Fig. ^ represents the experimental 
values of (69,61) to within one standard deviation inferred from pionic hydrogen data. The ultimate (60,61) values 
that are consistent with both the pionic hydrogen and the pionic deuterium data fill the area of the black strip. 



III. FINITE RANGE APPROACH 



Thus far our treatment of the pion-deuteron scattering problem has been carried out exclusively within the zero- 
range model. Although this model has served us well, it is based on certain idealization whose validity and conse- 
quences need to be examined. We therefore turn now to the question of formulating a finite range version of the 
approach presented in the preceding section. Relaxing the zero-range limitation has of course its quid pro quo in that 
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we have to worry now about the ofF-shell extension of the ttN scattering amphtude and this means that the pionic 
hydrogen problem has to be considered ab intitio in order to provide the necessary input for the 7rd calculation. 
Anticipating the application in the Faddeev type calculation, it will be convenient for us to work with separable 
potentials. To get insight into the pionic hydrogen problem, let us consider a two-channel situation, where the upper 
channel labeled as 1 corresponds to the neutral Tr^n system and the lower channel labeled as 2, respectively, to the 
■K^p system. We assume that the two-channel interaction respects isospin invariance and the isospin symmetry is 
broken only by the Coulomb potential operative in channel 2 and by the mass splitting within isospin multiplets. 
Since we wish to consider an atomic system it is essential to treat the Coulomb interaction exactly. To meet this re- 
quirement, we choose the two-channel Lippmann-Schwinger equation as our dynamical framework that in coordinate 
representation takes the form 

ui{r)^ / {r\G+{W)\r')[Vii{r',r")ui{r") + Vi2{r'y')u2{r")]dr' dr" (27a) 
Jo 



U2{r)^ / {r\G+{W)\r')[V2i{r',r")ui{r") + V22{r'y')u2{r")]dr' dr" (27b) 
Jo 

where we have assumed spherical symmetry of the problem and Uj (r) denotes zero orbital momentum radial wave 
function in channel j . The strong ttN interaction is adopted here in the form of a non-local potential matrix Vij . 
In (|2^) we have introduced the Green matrix whose only non-vanishing diagonal elements are 

(r|G+(W^)|r') =-(2A^l/pl)exp(^plr>) sin(pir<), (28) 

for the neutral channel, while in the charged channel we have to take into account the Coulomb interaction and the 
exact Green's function in this case reads 

{r\Gt{W)\r') = -{2^l2/p2)[Go{v,P2r>) + iF^{rj,P2r>)] Foiv,P2r<), (29) 

where r< = min(r, r'), r> ~ max(r, r'). In (p7|)-(p9|) W denotes the total cm. energy (including the rest mass), 
/ij are the reduced masses in the two channels and pj are the channel momenta: pj = zL^2fj,j{W — Ej) with Ej 
being the threshold energies and the sign ambiguity will be resolved in a moment. All masses here are assumed to 
take their physical values. In (p9[) r] ~ — a/i2/p2 and Go, Fq denote the standard Coulomb wave functions for orbital 
momentum £ = 0, defined in p6| . Finally, it should be noted that there is no ingoing wave in (p7|), as appropriate 
for a bound state problem. 

As mentioned above, to simplify matters, we assume that the interaction is separable, i.e. that the potential matrix 

is 

V,j{ry)=-v{r)s,,v{r'), (30) 

where the function v{r) represents the shape of the potential and the dimensionless parameters Sij are the measure 
of the strength of the potential. Time reversal implies Sij = Sji. With separable potentials, the system of integral 
equations (^^ can be solved analytically. To this end it is sufficient to multiply each of the equations by v(r) and 
integrate over r. This gives a system of two homogeneous algebraic equations for the two unknown quantities 

/•OC 

X, - / v{r)u,{r)Ar, j-1,2 
Jo 

and the latter will have a non-trivial solution if, and only if, the determinant of the system D{W) vanishes. Expanding 
the determinant, we are led to the explicit bound state condition 

D{W) - (1 + {v\Gt{W)\v)) (1 + S22 {v\Gt{W)\v)) - s?^ {v\Gt{W)\v) {v\Gt{W)\v) = (31) 

where we have introduced the abbreviation 

{v\GUW)\v) = / v{r) {r |G+(I^)| r') v{r') dr dr'. 
Jo 

The determinant can vanish only at some particular value of the energy W = Eb that will be interpreted as the 
bound state energy. Normally, knowing the underlying interaction, by solving ( pl| ) one obtains the binding energy. 
However, in the problem at issue we have a reversed situation: we know the binding energy from experiment and it 
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is the interaction that we are after. In the case of the pionic hydrogen atom we have an unstable bound state in the 
charged channel and the binding energy will be a complex number. We set 

Eb^E2+ - (e + zi7) (32) 

where Eis = — /i2a^/2 is the purely Coulombic Is state binding energy. Since in our formalism there is no room 
for the radiative decay of the pionic hydrogen the partial width 7 is a fraction of the total width F given by the 
formula 7 — r/(l + P^^) where P is the Panofsky ratio. It has been shown in ref. that the effect of the (7r~,7) 
reaction on the accounted for hadronic channels is negligible. The experimental values for e,r (cf. [Q) and P (cf. 
psf) adopted in this work, are 

e = 7.108 ± 0.013(stat) ± 0.034(syst) eV, 
r = 0.868 ± 0.040(stat) ± 0.038(syst) eV, 
P = 1.546 ±0.009, 



and in the following we shall take e = 7.108 ± 0.047 and 7 — 0.527 ± 0.047 as the input values. It must be 
immediately explained here that in this work we have defined e in accordance with a different convention, so that 
our e has opposite sign than that used in ref. [Q. In our approach we have tacitly assumed that under perturbative 
treatment all electromagnetic corrections contribute the same amount to the purely Coulombic level and to the level 
shifted by strong interaction. More precisely, we are going to ignore the small effects caused by the distortion of 
the wave function. Accordingly, the electromagnetic corrections need not concern us here and they have been left 
out altogether but, of course, they would be indispensable for calculating the total displacement of the level from its 
Coulombic position. 

The pole of the T-matrix that corresponds to the solution of ( pl| ) can be located on one of the four Riemann 
sheets as appropriate for a two-channel problem. This is also apparent from the mentioned above sign ambiguity in 
the definition of the channel momenta in (|2^). The right choice of the Riemann sheet is essential and this can be 
accomplished by proper adjustment of the signs of the imaginary parts of the channel momenta pj. We are using 
here the standard enumeration of the Riemann sheets introduced in ref. [E9l, i.e. 



sheet I 
sheet II: 
sheet III 
sheet IV 



Impi > 
Impi < 
Impi < 
Impi > 



Imp2 > 
Imp2 > 
Imp2 < 
Imp2 < 0. 



In the pionic hydrogen case, with an unstable bound state in channel 2, we have to enforce the pole to be located on 
the second sheet. 

To proceed further we need some concrete shape factor v{r) and our choice here is the exponential shape, i.e. we 
set 



;(r) = y/WJJj' exp(-/3r) 



(33) 



where fj, is the reduced pion-nucleon mass in the case of exact isospin symmetry (we take average mass for each isospin 
multiplet) and (3 is the inverse range parameter. With the exponential form (^3|), the potential (^) is identical with 
the familiar Yaniaguchi potential and the Green's function matrix elements can be obtained in an analytic form. The 
final result is 



{v\G+{W)\v) 



1 



fi {i-ipi/PY 

for the neutral channel, while the corresponding formula for the charged channel reads 



1 



{v\Gt{W)\v)=-^-^- 

n [i-iP2/py 



2Fi(l, iTj; IT) + 2- z^) 



1 



(34) 



(35) 



with z = (/3 + tp2)/{l3 — ip2)- The last fraction in (|3g) accounts for the Coulomb interaction and the symbol 
2-^1(0, 6; c^z) denotes the hypergeometric function defined in p^ . The computation of the hypergeometric function 
entering (35) is greatly simplified owing to the fact that the first parameter is equal to unity in which case the 
continued fraction representation of 2-Fi(l, c; z) discovered by Gauss proves to be useful. The continued 
fraction summation converges in the whole of the complex z plane away form the branch cut on the real axis running 
from one to infinity. 
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With exact isospin symmetry the three strength parameters sn, S12, S22 are not independent and can be expressed 
in terms of isospin 1/2 and isospin 3/2 strengths denoted hereafter as Si and S3, respectively. In the bound state 
condition both the real and the imaginary part of D{Eb) have to vanish simultaneously and that gives us two 
real equations. Since the bound state energy is known (cf. (p2|)), we put su = (si +2s3)/3 ; S22 = (2si + S3)/3 ; S12 = 
\/2(s3 — si)/3 in (|3l|), and regarding si and S3 as our two unknowns, we arrive at two algebraic equations that can 
be solved analytically 

Si Im ac* + si Im(o6* — c) — Im b — 0; (36a) 



S3 = -(l + siRea)/(Refe + siRec), (36b) 

where a = {{v\G+ (Eb^v) + 2 {v\G+iEB)\v))/3;b = (2 {v\G+ iEB)\v) + {v\G+ iEB)\v))/3 and c = 
(^v\Gi {Eb)\v') (^v\G2 {Eb)\v') . With si and S3 in hand, the corresponding scattering lengths {a2i with 1=1/2 and 
3/2) are obtained from 

a2i = -r- . (37 

P 1 - 52/ 

For local potentials the method outlined above could be also applied but in such case it would be more convenient 
to use instead of (|2^) an equivalent set of two coupled Schrodinger equations. For fixed energy and the proper choice 
of the Riemann sheet, these differential equations can be integrated numerically and the bound state equation is 
obtained from the requirement of vanishing of the Wronskian determinant. The latter is again a function of the 
isospin 1/2 and isospin 3/2 strength parameters, or if one prefers, the corresponding potential depths. Although the 
bound state condition is defined then only numerically but from it one can get two real equations that can be solved 
numerically using standard procedures. With a local potential, however, the solution of the three-body problem 
becomes much more complicated and this is the main reason why we preferred to work with a separable potential. 

Our calculational scheme is now complete and we shall present our results. Using as our input the experimental 
values of the pionic hydrogen level shift and width, the bound state equation has been solved analytically by adopting 
a number of "reasonable" values for f3 and in our computations we have used the values from 2 fm~^ to 10 fm~^. 
Although, we do not know the precise value of the range but there is no physical mechanism known that might 
generate long range forces in the ttN system, the longest range is unlikely to be bigger than 0.5 fm and this sets 
the lower limit of acceptable /3 values. In principle, there is no upper limit for /? but for /3 > 10 fm~^ we have in 
practice reached the limit of the zero-range forces and things change very little above that limit. The exact solutions 
of the bound state equation are presented in Table ^ where the errors reflect only the experimental uncertainty of 
our input, i.e. e,r and the Panofsky ratio. Our isoscalar and isovector scattering lengths are in good agreement 
with the values extracted in This has been illustrated in Fig. ^ where we have compared a representative sample 
of our computations with the values obtained by Sigg et al. The solutions corresponding to /3 spanning the 
range 2—10 fm~^ are located very close to each other in the (bi, bo) plane and putting more than three points 
on the plot might have obscured the picture. The error bars reflect only the experimental uncertainty of our input. 
As mentioned above, the bound state equation (pj) yields a second order equation for si and S3 and therefore we 



have always two solutions (cf. (|36D). Only one of them is presented in Table [V whereas the second solution leads 
to both 60 and 61 positive and has had to be rejected. When the two strength parameters si and S3 are known 
we can calculate not only the scattering lengths but also the effective ranges in each of the two isospin states and 
these values are presented in Table Instead of the effective range we use the parameter B21 that is defined from 
the expansion of the real part of the s-wave scattering amplitude in powers of the cm. momentum k, i.e. close to 



threshold, we have Ref2i{k) = 02/ -I- B21 k + ■ ■ ■ . For comparison, at the bottom of Table IV we give the values 



of all parameters inferred from a recent phase shift analysis The calculated scattering lengths, listed in Table 



[V, are almost independent upon (3, in contrast with the slope parameters B21 which change quite a bit when (3 is 
varied in the interval 2-10 fm~^. In addition to that, our values turn out to be always positive and therefore 
have opposite sign than those deduced from phase shift analysis . Actually, the pionic hydrogen data provide a 
strong constraint only for the scattering lengths and sticking to a simple ttN Yamaguchi potential it is not possible 
to get -B3 negative just by varying (3. Indeed, for fixed a2/ the slope parameter B2/ is given by the exact formula 

and since the expression in the square bracket is necessarily positive the sign of B2/ is bound to be opposite to that 
of a2/. To obtain a negative B3 a more sophisticated potential involving both repulsion and attraction would have 
been required There is no need for such extension, however, because our model has been devised for describing 
only the near threshold phenomena and is quite adequate at that. Expanding the phase shift close to threshold in 
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powers of k, we have S21 = a2ik + 0{k^) and it is apparent that a model providing merely the scattering length 
reproduces satisfactorily the phase shift in the neighbourhood of zero where 621 exhibits a linear behaviour. In our 
case this is all that counts as we never deal with higher energies. This means that the determination of the slope 
parameters is out of reach within our model since the appropriate energy scale has been set by the Coulomb energy 
in the pionic hydrogen, in which case terms proportional to B21 make negligible contribution. For an assigned value 
of /3 the slope parameters may be calculated but they are of no physical significance and comparing them with those 
resulting from phase shift analysis does not make much sense. 

As noted in 7|j27|], at the energy value close to the unstable bound state in channel 2, the scattering amplitude in 
the open channel 1, shows a strong resonant behaviour. For a separable potential, the s-wave scattering amplitude 
f{W) in channel 1 may be easily calculated analytically and takes a simple form 



fiW) = e smVp, = -- (^i+pyp2)2 ^(^) > (38) 



where S is the corresponding phase shift that for real W below the Tr^p threshold is a real number. The resonance 
is not of a Breit-Wigner shape but its position Er may be easily established from (^8|) as the energy at which the 
phase shift is equal to ^tt. Close to the resonant energy, i.e. for W ^ Er we have cot (5 « {W — Er)/{^Tr) and this 
allows us to infer the value of the width F^ of the resonance. In ref. Q the values of (e, 7) have been calculated by 
identifying them with {E2 + Eis — Er,Tr). In principle, the values of (e,7) obtained that way do not have to be 
identical with those determined from the position of the bound state pole. To check that point, we have repeated 
the procedure of ref. |7| but using our separable potentials whose depths have been adjusted to reproduce the values 
of (e,7) obtained in jl]. We found that the two methods give nearly identical results and the differences in (e,7) 
did not exceed 1 meV. For illustration, in Fig.^ we show the behaviour of sin S close to the resonance for the case of 
P — 3 fm~^ where the strengths parameters inferred from the pole location were si — 0.271820 and S3 = —0.245868. 

Before concluding our discussion of the pionic hydrogen we wish to mention one last thing, namely we are going 
to show how from (^l|) one can retrieve the Deser-Trueman formula (cf. This task will be accomplished by 

obtaining an approximate solution of (^) and to this end (31) is cast to the form 



l + Sos{W) {v\Gt{W)\v) ^0, (39) 
where we have introduced an effective energy dependent complex strength parameter ScS, defined as 

SosiW) = S22 - sj, {v\G+{W)\v) I (1 + sn {v\GX{W)\v)) . (40) 
The complex 7r~p scattering length a^rp can be expressed in terms of Scft(T^) evaluated at threshold 

_ Scff(S2) 



/1/3 1 - Scff(-E2) ' 



(41) 



and the Coulomb corrected tt p scattering length, denoted as a^^, can be obtained from the exact formula derived 
in Hi 

1/a^p = eVa^P + 2/i2a Ei(e), (42) 

where ^ — 4a^2/i5 and Ei(^) is the exponential integral function defined in |2^. It should be noted here that the 
zero-range limit (/3 ^ 00) does not exist in ( ^2| ) because the function Ei(^) for ^=0 has a logarithmic singularity. 
For the case of /3 = 3 fm~^ just considered, we obtain 

a^p = 0.12081 + 10.004441 fm; 
alp = 0.12068 + ^ 0.004458 /to; 

so that the Coulomb corrections do not exceed a fraction of a percent. However, in general, the Coulomb correction 
is model dependent, and, in particular, it is rather sensitive to the range of the nuclear potential what can be seen 
when the above result is juxtaposed with the 7rd case where the range of the potential was comparable with the size 
of the deuteron and, accordingly, the Coulomb correction to 7rd scattering length was much bigger (1.5%). 

Since we wish to obtain an approximate solution of (|3^) that is located not far from the Coulomb bound state, we 
set W = E2 + Els + SE where 6E is a small displacement. To calculate SE and derive the Deser-Trueman formula 
from (p9), we have to assume that (i) the complex energy shift SE = —e — «^7 is small in comparison with Coulomb 
eiieTgy(\SE / Eis\ « 1), and, (ii) that the range of the strong interaction is small as compared with the Bohr radius 
(/3 >> iJL20i). Introducing a complex momentum pc = y/2jI^^E\l = ifi2a corresponding to the Coulomb bound state, 
we can see that when p2 — > Pc then irj ^ —1 and the Green's function (pa) occurring in (p9h becomes singular. This 
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singularity is of paramount importance since it induces a zero in the nuclear S-matrix that is necessary to cancel the 
bound pole in the Coulomb S-matrix. As a result of this cancellation, the full S-matrix in the charged channel, which 
is a product of the Coulomb S-matrix and the nuclear S-matrix, remains finite at p2 = Pc- In compliance with the 
small shift assumption, we set p2 — Pc + ^P where 6p is supposed to be a small correction and since the most rapid 
variation in (js^) arises on account of the pole term, we approximate 1 + trj hy Sp/pc- Apart from that, elsewhere 
we replace p2 by Pc- The hypergeometric function for trj = —1 reduces to a polynomial 1 — and neglecting small 
terms of the order of pc/ /3, from (|3^) we obtain 

Sp w -Ai {pI/13) i^i2/^^)ses{E2) ~ -2ipl a^p (43) 

where we have used ( ^ ) retaining only linear term in a^rp- The above result gives the Deser-Trueman formula pd| ] 
in its standard form 

5E « Pc Sp/n2 « ot" a^p, (44) 

where, in view of the above discussion, it does not really matter whether we take a^p or a^p. It is perhaps in order 
to recall that although the Deser-Trueman formula (|4^ ) has been derived here for a specific choice of the underlying 
interaction, but its validity is quite general. To examine the accuracy of Deser-Trueman formula we turn again to 
our previous example when /? — 3 fm^^ and by computing a^p from ( ^l]) and inserting in (^), we obtain (e,7) = 
(7.024, 0.516) eV to be compared with our input values equal (e,7) = (7.108, 0.527) eV that ought to have been 
reproduced if formula (^) had been exact. It is a remarkable property of the Deser-Trueman formula that it is 
independent of the range of the underlying interaction and therefore the error in this formula must be of the same 
size as the uncertainty in the exact result caused by varying /3. If one is prepared to tolerate such uncertainty formula 
(4^) could be used to infer oi and 03. Introducing a two-channel K-matrix, isospin invariance can be invoked to pin 
down its elements at the single unsplit threshold 




and the complex tt p scattering length takes the form 

anp = K22 + iptKl2/{l -iptKii), (45) 



where pt is the momentum in the ir^n channel evaluated at the n^p threshold. The scattering length (45), unlike 
(^), does not depend upon the range. Inserting (^5|) in ( ^4[ ) and separating the real and the imaginary part, we end 
up with two real equations for the two unknowns ai and 03. To more than sufficient accuracy, the explicit solutions. 



x±y{l-2pty)/^j2^ /(l-pty)] (46a) 



03 = [a; T 2/(2 - Pty)/\/2iHy\ /(I - Pty), (46b) 

where x = e/2/X2a'^,?/ = ^j/2^2'^^ ^-nd the double sign in ( ^6| ) stems the fact that eq. ( p5| ) is quadratic in 02/. If 
(e,7) have been obtained in a model independent way then the results (^) are also model independent. As seen 
from Table ^ the uncertainty on ai and 03 (3% and 9%, respectively) induced by experimental errors on (6,7) is 
much bigger than the uncertainty caused by varying /? (about 1%). Under these circumstances it is perfectly justified 
to infer the ttN scattering lengths via Deser-Trueman formula and their numerical values obtained from (^6|) are 
displayed in Table tV whereas the corresponding 60 and 61 are presented in Fig. ^ 



It is apparent from Table IV that to improve upon Deser-Trueman formula we need some additional clue concerning 
P and it becomes something of a challenge to find ways to ferret out more precisely what the value of /3 might be. 
So far in our considerations we have not mentioned yet the pionic deuterium data and at this stage it is logical 
to ask whether this additional information might not help to pin down the range parameter of the ttN potential. 
Therefore, in the next step, we use the values given in Table as input for a three-body calculation, i.e. we use 
the separable potential (^3|) in the Faddeev equations for calculating the ird scattering length. The results of our 
computations are presented in Fig. ^ where we have plotted the ird scattering length versus p. The full circles 
represent the results obtained by the including the p-wave interaction (more precisely, just the P33 wave), while the 
open circles correspond to a situation where the delta has been left out. For reasons of clarity of the presentation 
these two sets of points have been given at different /3 values. The indicated error bars reflect the uncertainty in the 



input values (cf. Table IV ). For comparison, the experimental value of nd scattering length to within one standard 
deviation is given in Fig ^ as the area between the two horizontal lines. The striking feature apparent from Fig ^ 
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is that the resuhs are almost independent of the range parameter /3. Furthermore, the calculated scattering lengths 
are consistent with experiment for all /3 no matter whether the delta has been included or not. This result may come 
as a disappointment since the deuteron data give no illumination how to bracket the value of f3. 

In order to understand how the above result comes about we shall invoke again the static model, taking advantage 
of the fact that with the Yamaguchi potential representing the ttN interaction the static solution of the Faddeev 
equations may be readily obtained (cf. ref. ). Thus, introducing the Yamaguchi form factors and going to 
the static limit we can repeat the procedure outlined in the preceding section. The static solution of the Faddeev 
equations may be then sought in the form 



771 0^ 



and the above ansatz used in the Faddeev equations yields a set of two integral equations that differ from (|T^) in 
that the appropriate kernels contain now an extra factor 1/[1 + {k + k'Y / 0^]"^ . Despite this additional complication, 
the Fourier transform of this extended kernel still can be effected and leads to a simple analytic expression 



k'^ + {k + k'Y [P"^ + {k + k'YY J ^ T \ ' 2^ 13 

Using the above formula, similarly as before, we end up with a system of two algebraic equations for A(r) and X{r). 
Neglecting the binding energy correction (k — > 0), the resulting equations differ from (^) in that the zero-range pion 
propagator 1/r has to be multiplied by the function g[r) given by the formula 

,g(r) = l-e-'3'^(l + i/3r). (47) 

Therefore, the sought for solution for A{r) follows from ( |l0| ) after replacing 1/r by g{r)/r. Formula ( ^ ) proves to be 
quite useful for estimating the size of the /3 dependent correction and to this end we need to evaluate g{r) at some 
average value of r and a plausible candidate for such average value is the deuteron radius = V < r"^ > ~ 2 fm. 
Indeed, with this choice the second order formula ( [ll] ) that provides for a major contribution to OTrd will be little 
affected since by setting r = r^, we get (1/r) = 0.5 fm^^, not far from the values listed in Table ||. When /3 is varied 



in the range 2—10 fm^^, we have idP >4 in the exponential damping factor in (47), so that the P dependent terms 
make a contribution to g{r) at the level of a few percent and the resulting 7rd scattering length is almost independent 
upon p. This feature, sustained in the full Faddeev solution, is a consequence of the fact that the adopted range of 
the ttN forces was small as compared with the deuteron radius. 

In conclusion, we have seen that the uncertainty in the calculated ai and 03, as well as in OTrd, connected with 
the lack of knowledge of the range parameter constitutes only a small fraction of the uncertainty resulting from the 
experimental errors on the pionic hydrogen data. The above results may be viewed as an a posteriori justification of 
our zero-range model developed in Sec. O: introducing a finite range would be merely a fine tuning which is not yet 
affordable in the current state of affairs. 



IV. DISCUSSION 



Assuming that the underlying ttN interaction is isospin invariant, we have analysed the recent pionic hydrogen and 
pionic deuterium data with the purpose to extract from them ttN s-wave scattering lengths 02/ for 1=1/2 and 1=3/2. 
It is an empirical fact that the complex energy shift in either of these two atomic systems is small when compared 
with the corresponding Coulomb energy and with the appropriate Bohr radii setting the length scale, the 7r-p and 
TT-d interactions are of a short range. Under these circumstances Deser-Trueman formula provides an extremely 
good approximation, relating in a model independent way the Is level shifts and widths in the pionic hydrogen and 
pionic deuterium to the complex scattering lengths a,rp and OTrd, respectively. However, to infer 02/ from the latter 
quantities is a non-trivial dynamical problem and to be able to solve it we introduced a simple and transparent 
potential representation of the ttN interaction. Within this model we obtain explicit solution of the 7r~p bound state 
problem and also of the related three-body 7rd scattering problem at zero energy. 

We have assumed throughout this work that the ttN forces are of a very short range and this supposition follows 
from a particle exchange picture: there is no sufficiently light particle presently known that might be capable of 
generating forces whose range would exceed 0.3-0.4 fm (which roughly corresponds to a vector meson exchange). 
In this situation it was logical to take the zero-range limit as our point of departure. In order to find out what 
the deuteron data can teach us about ttN scattering lengths, we calculated the 7rd scattering length by solving the 
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appropriate three-body ttNN problem. This task was accomphshed, both within the static approximation, and also 
by using the Faddeev formalism. We demonstrated that the same static formula for a,rd can be derived from: (i) a 
set of boundary conditions; (ii) a static solution of Faddeev equations, and (in) a summation of Feynman diagrams. 
The static formula expressing a^rd in terms of ttN scattering lengths was found to be surprisingly accurate: the error, 
estimated by comparing the static result with the full Faddeev solution, was at the level of 2%, i.e. of the same size 
as the experimental error on a^rd- The standard second order formula was shown to be insufficient: the incurred error 
was three times bigger than the present experimental uncertainty on a,rd- Using as input the ttN scattering lengths, 
that had been inferred earlier ||7| from pionic hydrogen data, we obtained Ottc; by solving the Faddeev equations 
for zero-range ttN forces. The requirement that the calculated ajrc; be in agreement with experiment to within one 
standard deviation, imposes bounds on the isoscalar and isovector ttN scattering lengths. The values of the ttN 
scattering lengths determined that way, consistent with both the pionic hydrogen and the pionic deuterium data, are 
presented in Fig. |l]. 

In the next stage of this investigation we lifted the zero-range limitation introducing a range parameter. The pionic 
hydrogen bound state problem was solved afresh for a variety of range values. We derived the appropriate bound state 
condition and taking the Is level shift and width of the pionic hydrogen as input, we used this condition to determine 
the s-wave ttN potentials. This was possible since a complex condition is equivalent to two real equations, which 
for an assigned range, can be exactly solved for the 1=1/2 and 1=3/2 depth parameters entering the ttN potentials. 
Knowing the potentials, it was a trivial matter to calculate the corresponding s-wave scattering amplitudes. As can 
be seen from Table the resulting ttN scattering lengths are rather insensitive to the adopted value of the range 
parameter. 

The analysis of the pionic hydrogen presented in this work parallels that given in [0. We differ, however, in the 
adopted dynamical frameworks: in 0| Klein-Gordon equation together with a local ttN potential has been used, 
whereas we consider a non-relativistic Lippmann-Schwinger equation (equivalent to a Schrodinger equation) with a 
separable ttN potential. As may be seen from Fig. ||, the ttN scattering lengths inferred in this paper are in good 
agreement with those deduced in Q . This is a direct consequence of the fact that Deser-Trueman formula provides 
such a good approximation that we can make considerable progress in deducing the ttN scattering lengths without 
committing ourselves in great deal to the nature of the ttN dynamics. Since Deser-Trueman formula depends neither 
upon the shape of the ttN potential nor upon its range, the small changes in the ttN scattering lengths caused by 
varying the range parameter, must be attributed to the differences between the approximate Deser-Trueman formula 
and the exact range dependent solutions of the bound state equation. Thus, Fig. || illustrates the accuracy of 
Deser-Trueman formula. 

For an assigned range value, the pionic hydrogen data specify completely the ttN potentials, so that they may 
be used in the Faddeev equations in order to obtain the 7rd scattering length. The latter quantity was shown to 
be almost independent upon the range parameter (cf. Fig. ^) but was rather sensitive to the values of the ttN 
scattering lengths used as input in the Faddeev equations. The above finding, supporting the zero-range approach, 
could be explained by the fact that the range of the ttN interaction that was considered physically justified was small 
in comparison with the deuteron size. 

We conclude that the lack of knowledge of the range of the ttN interaction is responsible for some uncertainty in 
the deduced ttN scattering lengths but this uncertainty is rather small, at the level of 1%. The main source of error 
is still the experimental uncertainty in the pionic hydrogen data. 

It is rather obvious that the presented model contains several omissions but we think that they are not too severe, 
especially that the investigation has been confined to near threshold phenomena. As in all non-relativistic models 
based on static potentials virtual particle production, crossing symmetry, retardation and relativistic effects have not 
been even touched upon. Besides that, a separable potential is not considered to have a strong theoretical basis and 
has been adopted here merely for convenience as it simplifies considerably the solution of the Faddeev equations. 
There are also limitations on the completeness of the Faddeev approach where by restriction to three-body channels 
we were forced to leave out a wealth of inelastic features. The absorption channels leading to two-nucleon states 
are not easily incorporated in a Faddeev theory and require considerable enlargement of the present model which 
does not seem to be currently justified. While cognizant of the above deficiencies, we wish to believe that they are 
outweighted by the model merits. 
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TABLE I. 7rd scattering length obtained from the static model and from a Faddeev calculation in the zero-range model for 
different 60 and fei. For the NN forces we used PEST and Yamaguchi potentials, the results for the latter case are presented 
here in brackets. All entries are in 10~^/m^ units. 



61 


60 


model 


-9.47 


-9.05 


-8.63 



-0.65 



2-nd order 
static (IC) 
static (2C) 
Faddeev 
ditto with A 



-4.22 
-3.89 
-3.44 
-3.97 
-3.59 



■4.87) 
■4.21) 
■3.77) 
■4.27) 
■3.97) 



-3.97 
-3.69 
-3.29 
-3.76 
-3.37 



■4.57) 
■3.98) 
■3.58) 
■4.04) 
■3.73) 



-3.74 
-3.49 
-3.10 
-3.55 
-3.16 



■4.28) 
■3.77) 
■3.39) 
■3.81) 
■3.50) 



-0.22 



2-nd order 
static (IC) 
static (2C) 
Faddeev 
ditto with A 



-3.30 
-2.99 
-2.53 
-3.07 
-2.68 



■3.96) 
■3.32) 
■2.87) 
■3.37) 
■3.08) 



-3.06 
-2.78 
-2.36 
-2.85 
-2.46 



■3.66) 
■3.09) 
■2.68) 
■3.14) 
■2.85) 



-2.82 
-2.58 
-2.19 
-2.65 
-2.25 



■3.37) 
■2.87) 
■2.49) 
■2.92) 
■2.62) 



0.21 



2-nd order 
static (IC) 
static (2C) 
Faddeev 
ditto with A 



-2.38 
-2.08 
-2.62 
-2.16 
-1.76 



■3.04) 
■2.42) 
■1.97) 
■2.47) 
■2.20) 



-2.14 
-1.87 
-1.45 
-1.95 
-1.54 



■2.74) 
■2.19) 
■1.77) 
■2.24) 
■1.96) 



-1.90 
-1.68 
-1.28 
-1.74 
-1.34 



■2.45) 
■1.97) 
■1.59) 
■2.02) 
■1.73) 



TABLE IL The expectation values of r, 1/r and the values of vrd scattering length calculated for different NN wavefunctions. 
For ttN scattering lengths we have adopted their central values, i.e. bg = —0.22 and 61 — —9.05. All scattering lengths are 



given in 10 /m-n- units. 



NN wavefunction 





Hulthen 


PEST 


Paris 


Bonn 


(r) fm 


3.1345 


3.2309 


3.2685 


3.2536 


(1/r) fm-' 


0.55501 


0.45507 


0.44864 


0.46314 


2-nd order a^a 


-3.66 


-3.06 


-3.04 


-3.13 


static a^rd 


-3.09 


-2.78 


-2.78 


-2.82 


TABLE in. 7rd scattering lengths calculated from consecutive iterations of the Faddeev equations. 


All entries are in 


10 ^ /m^ units. 












PEST 


PEST 


Yamaguchi 


Yamaguchi 


order 


no A 


with A 


no A 


with A 


1 


-1.66 


-1.21 


-1.70 


-1.23 


2 


-2.98 


-2.66 


-3.42 


-3.30 


3 


-2.89 


-2.44 


-3.20 


-2.77 


4 


-2.85 


-2.48 


-3.11 


-2.91 


5 


-2.85 


-2.45 


-3.14 


-2.82 


6 




-2.46 


-3.15 


-2.87 


7 




-2.46 


-3.14 


-2.84 


8 






-3.14 


-2.86 


9 








-2.85 


10 








-2.85 
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TABLE IV. ttN scattering lengths inferred from pionic hydrogen data {B21 are slope parameters defined in the text). 



P 


a 1 


as 


Bi 


B3 


[tm J 




[iu m^ J 


[iU m^ J 


[iU m^ J 


z.u 


U.l/D/zt U.UU40 


-u.yo/ / =t u.uooz 


-D.DO 


i.yo 


3.0 


0.1760± 0.0045 


-0.9306± 0.0846 


-3.60 


0.81 


4.0 


0.1757± 0.0046 


-0.9263± 0.0841 


-2.46 


0.43 


5.0 


0.1756± 0.0046 


-0.9228± 0.0837 


-1.90 


0.27 


6.0 


0.1756± 0.0046 


-0.9197± 0.0834 


-1.57 


0.18 


7.0 


0.1756± 0.0046 


-0.9167± 0.0830 


-1.37 


0.14 


8.0 


0.1756± 0.0046 


-0.9138± 0.0827 


-1.23 


0.11 


9.0 


0.1756± 0.0047 


-0.9110± 0.0823 


-1.12 


0.09 


10.0 


0.1757± 0.0047 


-0.9082± 0.0820 


-1.05 


0.08 


Deser 


0.1760 ± 0.0046 


-0.9258 ±0.0857 






ref. i 


0.1679 ± 0.0059 


-0.785 ± 0.034 


-7.24 ± 3.06 


-4.08 ± 1.46 



FIG. 1. Constrains on the isoscalar and isovector scattering lengths imposed by pionic deuterium data. The black strip 
corresponds the one standard deviation region. The rectangle corresponds to the values obtained from pionic hydrogen data 
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FIG. 3. The values of isoscalar and isovector scattering lengths obtained by solving the bound state equation (^l|) for (3 
equal, respectvely, 2.0 fm~^, 6.0 fm~^, and 10.0 fm^^ (indicated on the plot). The point marked as Deser has been obtained 
from (^). The rectangle corresponds to the values obtained in [^. 
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FIG. 4. 7rd scattering length vs. the inverse range parameter /3 of the ttN potential. Full (open) circles correspond to a 
Faddeev calculation with (without) p-wave ttN interaction. 
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